Formation Histories of Dwarf Galaxies in the Local Group 
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ABSTRACT 

We compare the properties of dwarf galaxies in the Local Group with the simulated galaxies 
formed before reionization in a cosmological simulation of unprecedented spatial and mass reso- 
lution, including radiative feedback effects. We find that a subset of the Local Group dwarfs are 
remarkably similar to the simulated dwarf galaxies in all their properties already before reion- 
ization. Based on this similarity, we propose the hypothesis that Local Group dwarfs form in a 
variety of ways: some of them are "true fossils" of the pre-reionization era, some of them form 
most of their stars later, after reionization (we call them "survivors" of the reionization era), 
and the rest of them form an intermediate group of "polluted fossils" . We also identify a simple 
observational test that is able to falsify our hypothesis. 

Subject headings: cosmology: theory - galaxies: formation - stars: formation 



1. Introduction 

Nothing emphasizes our lack of understanding 
of the origin of the dwarf galaxies of the Local 
Group more than the well known "substructure 
crisis" of the Cold Dark Matter (CDM) paradigm 
(Klypin et al. 1999; Moore et al. 1999): the fact 
that the number of dark matter halos in the Lo- 
cal Group predicted by CDM simulations is much 
larger than the number of observed Galactic satel- 
lites. Several explanations have been suggested so 
far for resolving the "crisis" . 

The hypothesis of Bullock, Kravtsov, & Wein- 
berg (2000) - we will call it hereafter the "reion- 
ization scenario" - proposed that dwarf spheroidal 
(dSph) galaxies in the Local Group formed during 
the pre-reionization era and processes of photo- 
ionization feedback (Babul & Rees 1992; Efs- 
tathiou 1992; Shapiro, Giroux, & Babul 1994; 
Haiman, Rees, & Loeb 1996; Thoul & Weinberg 
1996; Quinn, Katz, & Efstathiou 1996; Weinberg, 
Hernquist, & Katz 1997; Navarro & Steinmetz 
1997; Gnedin 2000b; Dijkstra et al. 2004; Shapiro, 
Iliev, & Raga 2004, Susa & Umemura 2004) re- 



^ Institute of Astronomy, Madingley Road, Cambridge 
CMS OHA, England; ricotti@ast.cam.ac.uk 

^CASA, University of Colorado, Boulder, CO 80309, 
USA; gnedin@casa.colorado.edu 



suited in the suppression of star formation below 
the observable level in 90% of the low mass dark 
matter halos populating the Local Group. The 
"reionization scenario" is prompted by the fact 
that almost all Local Group dwarf Spheroidals ex- 
hibit a prominent old population and a decline 
of their star formation rates about lOGyr ago. 
While in the standard A-CDM theory (Spergel 
et al. 2003, Tegmark et al. 2004) stellar reioniza- 
tion took place 12.5 Gyr ago, current observational 
techniques for measuring absolute stellar ages are 
not able to differentiate between 10 and 12.5 Gyr 
(e.g., Krauss & Chaboyer 2003). Although the 
observed drop of star formation rate is consistent 
with reionization at redshift z 6, as noted by 
Grebel & Gallagher 2004, the metallicity spreads 
observed in the old populations may imply a star 
formation history protracted over about 2 Gyr, to 
redshift z ^ 3. In this paper we critically discuss 
this important issue. We will show that our simu- 
lations reproduce the observed metallicity spreads 
without a protracted period of star formation. 
The metallicity spreads observed in our simulated 
dwarfs can be understood in terms of hierarchical 
accretion of sub-halos containing stars with differ- 
ent metallicities. 

An alternative explanation to the "substruc- 
ture crisis" was proposed by Stoehr et al. (2002). 
According to their hypothesis, dwarf galaxies are 
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hosted in more massive dark matter halos 10^ — 
10^° Mq) than in the reionization scenario. While 
this idea resolves the "crisis" , it does not provide 
an explanation for the origin of dwarf galaxies. 

A more comprehensive and complex picture has 
been recently proposed by Kravtsov, Gnedin, & 
Klypin (2004). Under their hypothesis - which we 
will hereafter call the "tidal scenario" - most dwarf 
galaxies used to be much (up to a factor of 100) 
more massive during the formation episodes of the 
Milky Way and Andromeda, and so were not af- 
fected by photo-ionization feedback. They formed 
most of their stars at z ~ 3 but later most of the 
dark matter was tidally stripped during their evo- 
lution within the Milky Way and Andromeda ha- 
los. Kravtsov et al. (2004) emphasize that in their 
model a large variety of star formation histories 
is observed, and some of dwarf galaxies do indeed 
follow the "reionization track" , but the bulk of 
stars in the Local Group dwarfs are formed after 
reionization. 

Thus, the key difference between the reioniza- 
tion and tidal hypotheses is whether most of the 
stars in a given dwarf galaxy formed in halos 
of mass ;S 10® Mq before reionization {z > 6, 
12.5 — 13Gyr ago) or in more massive halos after 
reionization, during the epoch of formation of the 
Milky Way and Andromeda (i;'~2-3, 10-11 Gyr 
ago). 

In this paper we propose that dSphs are fos- 
sils of galaxies with mass Mdm < 10® M0 that 
formed before reionization while dirr formed later 
in more massive halos. The original contribution 
of our study is that for the first time we have been 
able to simulate the first population of galaxies 
including the relevant feedback processes that are 
known to self-limit their formation. Before our 
work (Ricotti, Gnedin, & Shull 2002a, 2002b) the 
mainstream wisdom was that negative feedback ef- 
fects (namely the rapid photo-dissociation of H2) 
suppress the formation of all galaxies with mass 
smaller than 10** solar masses. We find instead 
that these galaxies do form and, at the epoch of 
their formation, they closely resemble dSph galax- 
ies. This second result is the one we emphasize 
in the present paper. We also show that reioniza- 
tion feedback is not the dominant mechanism that 
suppresses star formation in most small mass ha- 
los. Negative feedbacks such as photo-heating by 
the stars inside each galaxy and inefficient cool- 



ing produce the observed properties of dSphs well 
before reionization of the IGM is complete. 

It is quite surprising that we find such an as- 
tonishing agreement between all properties of sim- 
ulated dwarf galaxies at 2: 8 and dSphs. Our 
simulations did improve on previous ones but are 
still far from a complete treatment of the problem 
that, as we will briefly illustrate, is non-trivial. 
The efficiency of star formation in the first mini- 
halos with masses ^10® Mq {i.e., virial temper- 
ature Tvir ^ 10'' K) is determined by many com- 
peting feedback processes, including the formation 
and disruption of H2, metal enrichment, photoe- 
vaporation and SN explosions. Since small mass 
proto-galaxies cannot cool by Lyman-alpha emis- 
sion, the presence of H2 molecules is essential to 
trigger their initial collapse. We know that H2 
is easily destroyed by FUV radiation emitted by 
the first stars (negative feedback). On the other 
hand ionizing radiation and dust production are 
known to enhance the formation rate of H2 (posi- 
tive feedbacks). Whether or not c£fic;icnt star for- 
mation is possible in these mini-halos depends on 
which feedback (positive or negative) is dominant. 
By including 3D radiative transfer of ionizing ra- 
diation (Ricotti et al 2002a) in the simulations we 
are able to model two additional feedbacks: the 
positive feedback for H2 formation and the photo- 
evaporation of dwarfs from internal sources. As al- 
ready mentioned, contrary to previous studies we 
find that low surface brightness galaxies do form 
in a fraction of the first mini-halos. Reionization 
does not affect substantially their properties be- 
cause they have already lost most of their gas due 
to internal ionizing sources. After reionization we 
expect that the formation of small mass galaxies 
will stop and their stellar populations will evolve 
only passively. Only the few most massive dwarfs 
in our simulation retained some of their gas and 
may continue to form stars after reionization. 

In summary the goal of this paper is to pro- 
vide evidences in support of the primordial origin 
of dSphs by comparing the abundant observational 
data on Local Group dwarfs with the results of de- 
tailed numerical simulation of formation of dwarf 
galaxies in the pre-reionization era. Aside for the 
origin of dSph galaxies the aim of study is to im- 
prove our understanding on theories of the forma- 
tion of the first galaxies and answer a longstanding 
question: how massive were the first galaxies? 
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2. Simulation 

The results shown in the present work arc based 
on a detailed analysis of the galaxies found in high- 
resolution cosmological simulations of a cosmologi- 
cal volume during the pre-ionization era. The sim- 
ulations include dark matter, gas, stars and the 
feedback of star formation on the gas chemistry 
and cooling. The main addition with respect to 
other works {e.g., Machacek et al. 2001, Tassis 
et al. 2003) is a self-consistent treatment of the 
spatially-inhomogcncous and time dependent ra- 
diative transfer of photons in four frequency bands 
(H lionizing radiation. He I and He II ionizing ra- 
diation and optically thin radiation). This allows 
us to study in addition to the global feedback ef- 
fects of the background radiation, also the role of 
internal and local feedback effects. Partially be- 
cause of the inclusion of a more realistic treat- 
ment of radiative feedback, we find that the H2 
photodissociating backgromid is not the dominant 
feedback and does not stop star formation in mini- 
halos. 

The simulation we use in this paper has been 

thoroughly described in Ricotti et al. (2002a, 
2002b) as run "256Llp3". Here we just remind 
the reader that the simulation included 256^ dark 
matter particles, an equal number of baryonic 
cells, and more than 700,000 stellar particles. 
The mass resolution of our simulation is 900 solar 
masses in baryons, and real comoving spatial res- 
olution (twice the Plummer softening length) is 
150/i~^pc (which corresponds to a physical scale 
of 24 parsecs at z = 8.3). This unprecedented 
resolution allows us to resolve the cores of all 
simulated galaxies that would correspond to the 
observed Local Group dwarfs. The simulation 
includes most of the relevant physics, including 
time-dependent spatially-variable radiative trans- 
fer, detailed radiative transfer in Lyman- Werner 
bands, non-equilibrium ionization balance, etc. 

In this work, however, we also include the ef- 
fect of rcionization in that simulation. Because 
the size of the simulation box has been fixed at 
lh~^ comoving Mpc, the simulation volume is too 
small to model the process of cosmological rcion- 
ization with sufficient accuracy. We, therefore, as- 
sume that the simulation volume is located inside 
an H II region of a bright galaxy at a higher red- 
shift. Specifically, we introduce a source of ioniz- 
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Fig. 1. — The mean gas fraction (thin lines) and 
the mean stellar fraction (thick lines) as a function 
of the dark matter mass for the simulated objects 
at 2; = 10.2 (dashed gray lines) and at the final 
moment of our simulation, z = 8.3 (solid black 
lines). Both fractions are normalized by the uni- 
versal baryon fraction. 

ing radiation within the computational box, prop- 
erly biased, which corresponds to a star-forming 
galaxy with the constant star formation rate of 
1 solar mass per year (similar to star formation 
rates of observed Lyman Break Galaxies at 2; ~ 4, 
Steidel et al. 1999). The source is switched on 
at z = 9.0, and by z = 8.3 the whole simulation 
box is completely ionized, and star formation in 
dwarf galaxies has ceased. Because of the compu- 
tational expense, we have not continued our sim- 
ulation beyond z = 8.3. As mentioned in the in- 
troduction the photoevaporation of dwarf galax- 
ies after rcionization has been studied extensively. 
Recent works include those of Shapiro, Iliev, & 
Raga (2004) and Susa & Umemura (2004) that ad- 
dressed the importance of self-shielding using 3D 
simulations coupled to ID radiative transfer {i.e., 
they assume that the ionizing source is very far 
from the dwarf galaxy). 

Figure 1 shows the mean gas (thin lines) and 
stellar (thick lines) fractions as a function of the 
dark matter mass for the objects in our simula- 
tion. We show these quantities at redshift z = 10.2 
(dashed lines) and z = 8.3 (solid lines). Notice 
that the objects with masses below about 10^ Mq 
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Fig. 2. — The stellar surface (filled circles) and 
volume (gray lines) densities (top row) and stel- 
lar (black lines) and dark matter (gray triangles) 
velocity dispersions (bottom row) as a function of 
radius for three representative objects from our 
simulation at z = 8.3. The degree of anisotropy 
of the velocity dispersion is usually measured by 
the parameter (3=1 — {(jr/(Jt)~'^ , where (T^ is the 
radial velocity dispersion and CTj is tangential ve- 
locity dispersions. The dashed lines in the bottom 
panels show Urjof 

contain no gas at the end of our simulation. Most 
of the small mass galaxies loose all their gas well 
before reionization that in our simulation happens 
between redshift z 9 and 8.3. This is the re- 
sult of gas photo-heating by massive stars inside 
the galaxy and not by the external ionizing back- 
ground. A remarkable result is that the mean 
stellar fraction is a time-independent quantity. It 
is apparent from the figure that the mean stellar 
fraction as a function of halo mass is well approx- 
imated by power law that does not evolve with 
redshift (the increase in the stellar fraction for the 
lowest mass objects is probably an artifact of the 
limited mass resolution). 

Figure 2 shows radial profiles for three represen- 
tative and well resolved objects. In the top panels 
we show surface (filled circles) and volume (gray 
lines) luminosity densities in the V band. The bot- 



tom panels show the velocity dispersions for both 
the dark matter (about 50, 000 particles per halo) 
and stars (about 10,000 particles per halo). No- 
tice that stellar velocity dispersions (black points) 
are below that of the dark matter (gray triangles), 
and that stellar orbits are mildly anisotropic in the 
outer parts. The dashed lines in the bottom pan- 
els show Or I Of, where Or and at are the radial and 
tangential velocity dispersions, respectively. 

Clearly, it is a long way to go from z = 8.3 to 
the present day. Therefore, in the rest of this pa- 
per, we are assuming that the stellar component of 
om- galaxies does not change after z = 8.3 except 
due to passive stellar evolution. This hypothesis, 
obviously, does not apply to all observed dwarfs, 
but it may apply to some of them - as we discuss 
later. The probability, P, of a galaxy collapsed at 
redshift 2;coii to survive to redshift z without being 
incorporated into a larger object is approximately 
P ~ (1 -I- z)/{l + 2coii) (Sasaki 1994). We can 
therefore estimate that about 10% of the galax- 
ies present in our simulation at z = 8.3 may have 
survived almost unspoiled until today. With this 
approach, we cannot determine what dark matter 
mass these galaxies would have by z = 0: they 
can continue accreting^ the dark matter and then 
lose a significant fraction of it during their orbital 
evolution within the parent galaxy halo (in this 
case - Milky Way or Andromeda). It is therefore 
important to keep in mind throughout this paper 
that we only use the properties of the stellar com- 
ponents of simulated galaxies. Thus, we cannot 
make a direct test of the tidal scenario, and have 
to rely on stars to tell us the story of the dark 
matter. 

In order to predict the luminosities and stellar 
densities of simulated galaxies at z = 0, we need to 
take into account stellar mass-to-light ratios and 
stellar evolution. Since the simulated stars would 
be old by z = 0, we adopt a value for the mass- 
to-light ratio Y as measured in old galactic glob- 
ular clusters, F ~ 1.5 — 2.0 (e.g., Mandushev et 
al. 1991, Pryor & Meylan 1993). In addition, after 

Note that a dark matter halo that virializes at z = 10 by 
the time it reaches redshift 2 = is about ten times larger 
in size and in mass even if it did evolve passively without 
accreting any satellite (i.e., in isolation). This is a conse- 
quence of the cosmological evolution of the Hubble constant 
and mean density of the IGM in which it is embedded (see 
Gunn & Gott 1972; Ricotti & Wilkinson 2004). 
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some 13 Gyr of stellar evolution, assuming Kroupa 
IMF we estimated (using Starburst99, Leitherer 
et al. 1999) that our stars will loose about 50% of 
their mass. The exact value depends on the as- 
sumed lower mass cutoff of the IMF. Thus, in or- 
der to compute the luminosities of the simulated 
galaxies at z = 0, we adopt a fiducial ratio of 
M^,{z = 8.3)/Lv{z = 0) = 3. However, our re- 
sults are not immensely sensitive to that number 
- changing it to 2 or 10 would not change our con- 
clusions (see Fig. 3 and Fig. 4). 

3. Results 

The core of the rcionization scenario is that 
some of the Local Group dwarfs are true fossils of 
the pre-reionization era. Clearly, not all dwarfs fall 
into this category, since many of them had recent 
episodes of star formation. But is there a subset 
of Local Group dwarfs that formed before reion- 
ization? We can attempt to answer this question 
by comparing our simulated dwarfs as they would 
appear at 2: = with the real observed galaxies. 
Figure 3 presents such a comparison for the ba- 
sic structural parameters: the luminosity Ly, the 
core radius Tc, and the central luminosity density 
Ic, defined as Ic = Ey/2rc, where Sy is the sur- 
face brightness. From the top panel of that figure 
it is clear that the simulated galaxies, shown by 
the gray squares without error-bars, do not exceed 
about 1O''L0. Moreover, it appears that the ob- 
served galaxies (shown by squares with error-bars, 
open circles, and filled circles - the different sym- 
bols will be explained later in the text) also fall 
into two groups - those with luminosities on the 
order of IQ^Lq and below, and those with lumi- 
nosities in excess of 5 x IO^Lq. In addition, these 
two groups are also clearly separated in Vc — Ic 
plane. The low luminosity dwarfs appear to trace 
the same distribution as the simulated galaxies, 
having smaller core radii as the central luminosity 
density increases. This observation can be under- 
stood in terms of the less efficient cooling of the 
gas in dark matter halos of smaller mass. 

Since our simulation has a finite size of the com- 
putational box, it is clear that we miss most lu- 
minous galaxies forming at the same time as our 
simulated dwarfs. The tilted dotted line roughly 
traces the upper envelope of the simulated galax- 
ies, so brighter galaxies that trace the same distri- 
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Fig. 3. — The V-band luminosity (top) and the 
core radius (bottom) vs the central luminosity 
density for the Local Group dwarfs (squares with 
error-bars, open circles, and filled circles - the dif- 
ferent symbols are explained in the text) and for 
the simulated galaxies passively evolved to z = 
(gray squares). Two dotted lines show two pos- 
sible boundaries between the "fossils" and "sur- 
vivors". The data is from Mateo (1998). We have 
complemented the data using updated distances 
and photometry from van den Berg (2000) and 
metallicities and metallicity spreads from Grebel 
et al. (2003). 



bution in the Ic — Ly plane would fall below that 
line. Some of the observed dwarfs are indeed in 
this region, but, as the bottom panel of Figure 3 
shows, their core radii are somewhat larger than 
those of the simulated galaxies. Thus, we can be 
confident that we do not miss a large fraction of 
more luminous dwarfs in our simulation that fol- 
low the same distribution as the faint ones. It does 
not, of course, mean that there should be no more 
luminous galaxies in a simulation with a larger 
volume - they will be present, but most of them 
will have masses in excess of 10^ M©. Therefore, 
because they will collapse via atomic line cooling, 
they should have a very different internal struc- 
ture and be much less susceptible to losing their 
gas after reionization. Figure 4 is the same as Fig- 
ure 3 but we assume M/Ly = 10 instead of the 
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Fig. 4. — Same as in Fig. 3 but assuming MjLy = 
10 instead of the fiducial value MjLy = 3. 
The agreement between simulations and data is 
mildly sensitive to the assumes mass to light ra- 
tio for all the scaling relationships except the 
luminosity-metallicity relation that gets worst for 
larger MjLy. 

fiducial value MjLy = 3. The purpose of this 
figure is to show that our results are robust and 
almost independent on the assumed M/Ly ratio. 

It is tempting, therefore, to hypothesize that 
the two groups of low and high luminosity dwarfs 
correspond to two distinct formation tracks - the 
low luminosity ones are the "fossils" of the pre- 
reionization era, and high luminosity ones formed 
their stars via atomic line cooling perhaps via 
Kravtsov et al. (2004) scenario, as tidally stripped 
more massive halos (in the absence of a better 
name, we will call them "survivors" of reioniza- 
tion era). 

We can now see how our separation in these 

two groups (based only on the luminosity) com- 
pare with the available data on the star forma- 
tion history of observed dwarfs. Hereafter we 
will adopt the star formation histories compiled 
by Grebel & Gallagher (2004). It is well known 
that almost all dwarfs in addition to a more or 
less prominent old population clearly show some 
more recent star formation. Therefore taken at 
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Fig. 5. — A tentative separation of all Local Group 
dwarfs as well as other nearby dwarfs into three 
groups: "survivors", that formed via the tidal sce- 
nario (mostly dirr), "true fossils" , that formed via 
the reionization scenario, and "polluted fossils", 
that went through both stages. In our tentative 
classification most dSphs (14) are true fossils: five 
are M31 satellites, six are Milky Way satellites 
and three (the underlined names) are isolated from 
both M31 and the Milky Way. 

face value, our hypothesis does not hold in com- 
parison with the data. However, it is unreasonable 
to assume that all of the galaxies that formed dur- 
ing the pre-reionization era formed absolutely no 
stars afterwards. Some of them should surely con- 
tinue accreting dark matter, becoming more mas- 
sive and switching from the "reionization" track to 
the "tidal" track. To be more specific, the ability 
of a dark matter halo to accrete fresh intergalac- 
tic gas is controlled by the filtering mass (Gnedin 
2000b). In the standard CDM-^A theory the fil- 
tering mass increases with time, as does the mass 
of an individual dark matter halo. Thus, the ques- 
tion of whether a given halo is able to accrete fresh 
gas after reionization (and, presumably, form new 
stars) depends on whether its mass grows faster or 
slower than the filtering mass. We can, therefore, 
expect that some of galaxies that formed before 
reionization will form no new stars, and some of 
them will form more stars as they fall into the 
Milky Way or Andromeda halo and their local fil- 
tering mass and accretion history vary. The spe- 
cific rate with which the gas is re-accreated must 
depend on the details of a satellite orbit, as has 
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been emphasized by Kravtsov et al. (2004) , and so 
a large variety of possible star formation histories 
is expected, in a qualitative agreement with the 
observational data. A much more advanced simu- 
lation than the one presented here will be required, 
however, to check whether a quantitative agree- 
ment between the theory and the data holds - for 
example, whether the relative fractions of dwarfs 
with different star formation histories are consis- 
tent with the data. 

As mentioned in the introduction another prob- 
lem faced by the "fossil" hypothesis is the large 
metallicity spreads observed in dwarfs with a dom- 
inant old population. Grebcl & Gallagher (2004), 
have suggested that these large spreads are not 
compatible with star formation protracted over a 
period of less than 1 Gyr (the time elapsed be- 
fore reionization) . They conclude therefore that 
star formation in old dwarfs was not affected by 
reionization. In defense of the "fossil" hypothe- 
sis we will show as part of our analysis that the 
metallicity spreads in our simulated dwarfs are 
consistent with observations. Note that star for- 
mation in these dwarfs was protracted for less than 
1 Gyr. Our result is probably due to accretion of 
stars from satellites with different metallicities and 
wouldn't be captured by any standard chemical 
evolutionary model. 

In order to schematically separate the observed 
dwarfs in groups according to their star forma- 
tion histories, we call those dwarfs that formed 
few (say, less than 30%) stars after reionization 
"true fossils" , and those that did form a sub- 
stantial amount of stars "polluted fossils". The 
various symbols in Figure 3 and the following 
ones correspond to these classes: "true fossils" 
arc shown with filled squares with error bars, 
"polluted fossils" arc shown with open circles. 
Galaxies with dominant younger stellar popula- 
tions (they mostly are dirrs) are shown with larger 
filled circles and we call them "survivors" . For the 
last two groups the error bars are not shown for the 
sake of clarity. All the galaxies in our simulation 
have almost identical properties to the "true fos- 
sils" and most "polluted fossils" . All these galax- 
ies have a prominent old stellar population as ex- 
pected. Instead all the galaxies with a prominent 
young population (filled circles) are part of a dis- 
tinct population that is not present in our simula- 
tion. 
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Fig. 6. — The V-band luminosity per cube of the 
core radius vs the central luminosity density for 
the Local Group dwarfs and for the simulated 
galaxies (gray squares) . In this and further figures 
filled black squares mark "true fossils", open cir- 
cles mark "polluted fossils" , and filled circles show 
"survivors" . The dotted line shows the power-law 
fit to the "true fossils" with the slope fixed to 3/2 
(see Gnedin 2000a for details). 

Figure 5 presents a separation of the Local 
Group dwarfs into these three types. Galaxies 
identified as "survivors" can be thought of as 
those dwarfs which formed after reionization (or 
mostly after reionization), and for which the to- 
tal mass exceeded the local filtering mass at some 
moment during their evolution. Polluted fossils 
formed a significant fraction of their stars in the 
pre-reionization era, but also exceeded the local 
filtering mass at some point in their history. The 
"true fossils" formed most of their stars in the pre- 
reionization era. As a weak argument in favor of 
this hypothesis, we notice that "polluted fossils" 
are on average brighter than the "true" ones. 

This separation of galaxies into individual 
classes is in some cases rather artificial. For ex- 
ample both Sextans and Fornax had some recent 
star formation. We labeled Sextans a "true fossil" 
because this recent episode was very small, while 
Fornax had a factor of 3 to 4 more stars formed 
after reionization, and so is labeled a "polluted 
fossil" . Based on their star formation history Leo 
A and Leo I could also be classified as low lumi- 
nosity "survivors". 

In order to support our conclusion, we consider 
other properties of the simulated and real galaxies. 
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Figure 6 shows the comparison between the data 
and the simulation for the relation between Ly/r^ 
and Ic, discussed in Gnedin (2000a). We find a 
scaling relationship 

shown by the dotted line. If the interpretation 
of Gnedin (2000a) of this relation is correct, then 
it simply reflects the Schmidt law for star forma- 
tion on scales of hundreds of parsecs - in that case 
the agreement between the data and the simula- 
tion is by construction, rather than a success of 
the model: the star formation in the simulation is 
assumed to follow the Schmidt law (1.5 slope in 
the relation between the star formation rate and 
gas density), and the amplitude is fixed by ob- 
servational constraints. On the other hand, this 
agreement does support the rcionization scenario, 
because "survivors" do not follow this tight rela- 
tion, which could be expected if their stars formed 
in several independent bursts. 

Note that L^/{Icr^) is the ratio of the total lu- 
minosity to the core luminosity and is therefore 
a measure of how concentrated the light of the 
dwarf galaxy is. This ratio is less than 10 for 
dSph galaxies like Draco and Ursa Minor {Ic ~ 
5 X lO~^i0pc~^) and is almost 100 for dEs like 
NGC205 and NGC147 {h ^ ILopc-^). The 
galaxies that we classified as "survivors" have gen- 
erally central luminosity density comparable to 
dEs but less concentrated light profiles. We do 
not observe such objects in our simulation - despite 
the fact that cores of these galaxies, if they were 
present in the simulation, would be completely re- 
solved. 

Figure 7 compares the surface brightness and 
the core radii of simulated and observed galax- 
ies as a function of the luminosity. The simu- 
lation reproduces very well the structural prop- 
erties and scaling relations of the galaxies that 
we classified as "true fossils". The lines show, 
for comparison, the scaling relations derived by 
Kormendy & Freeman (2004) for more luminous 
late-type galaxies (Sc-Im) with luminosities in the 
range lO^L© < is < lO^^L©. 

Figure 8 compares the dynamical state of these 
galaxies: the stellar velocity dispersion and the 
total mass-to-light ratio. The latter quantity for 
the simulated galaxies is computed in the same 
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Fig. 7. — Surface brightness (top) and core ra- 
dius (bottom) as a function of the V-band lu- 
minosity for the simulated and observed dwarfs. 
Symbol markings are the same as in Fig. 6. For 
comparisons we show (solid lines) the scaling re- 
lationships for the more luminous Sc-Im galaxies 
(IO^Lq <Lb^ 10"^©) derived by Kormendy & 
Freeman 2004. 

way as for the observed ones (M/L = 136a-^rc/ L, 
see Mateo 1998), and is not taken from the sim- 
ulations, because, as we have mentioned above, 
there is no simple way to relate the dark matter 
mass of simulated galaxies at z = 8.3 with their 
dark matter masses a.t z = 0. For a comparison, 
the solid lines show the same scaling relations for 
late- type galaxies (Kormendy & Freeman 2004). 
Within the errors the a — Ly relation seems to be 
the same power law for all the galaxy types, but 
"true fossils" obey a steeper relationship between 
the mass-to-light ratio and luminosity than dirr 
and Sc-lm galaxies. 

Figure 9 gives yet another comparison between 
the simulation and the data: the metallicity as a 
function of stellar luminosity. The agreement be- 
tween stellar metallicities in the simulation and 
the data is somewhat worse: while the spread in 
the distribution is consistent with the data, sim- 
ulated galaxies are more metal-rich by about a 
factor of 1.5 — 2. At this point we are not wor- 
ried by this disagreement as our simulation is not 
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Fig. 8. — Total mass-to-light ratio (top) and stel- 
lar velocity dispersion (bottom) for the simulated 
and observed dwarfs. Symbol markings and the 
meaning of the lines are the same as in Fig. 6. 
Fewer data points appear in this figure because 
the data do not exist for all Local Group dwarfs. 

sophisticated enough to predict the stellar metal- 
licities to within a factor of two. There are sev- 
eral physical and numerical effects that will load 
to overestimated stellar metallicities in the simu- 
lations, such as preferential ejection by supernovae 
of motal-cnrichcd gas (MacLow & Fcrrara 1999), 
imperfect mixing of the supernova ejecta on scales 
not resolved by a simulation, reduced gas accretion 
on the simulated galaxies due to a limited size of 
the computational box, etc. 

Low metallicity stars can be produced in two 
different ways: either by inefficient star formation 
combined with continuous gas loss, or in a single 
short burst. The first scenario is consistent with 
the properties of dSph galaxies which have an ex- 
tended low surface brightness stellar component. 
The second scenario is usually associated with a 
higher efficiency of star formation, which results 
in a more concentrated light profile than the one 
observed in dSph, and is probably more appropri- 
ate to explain the metal enrichment of globular 
cluster and dwarf Elliptical galaxies. 

We have already anticipated that our Simula- 
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Fig. 9. — The stellar metallicity for the simulated 
and observed galaxies. Symbol markings are the 
same as in Fig. 6. We adopt solar abundances and 
yield y=pz/p* = 1 Z©. 

tion reproduce the large metallicity spreads ob- 
served in dSphs. In Figure 10 we show a compar- 
ison between the simulated and the observational 
data for the metallicity spreads, (T[pg/^] , both as a 
function of the mean metallicity [Fe/H] (left panel) 
and function of luminosity Ly (right panel). 
We find that the typical metallicity variance in 
each dwarf is about 0.5 dex, in good agreement 
with observations. The star formation history in 
simulated dwarfs is characterized by multiple in- 
stantaneous starbursts extended over a period of 
0.5-1 Gyr. The starbursts often do occur almost 
the same time but in distinct dark matter sub- 
halos that are in the process of merging to form a 
more massive dwarf. The accretion of stars from 
the hierarchical distribution of sub-halos partially 
explains the large metallicity spread that we find 
in our simulation. 

Thus we conclude that the large metallicity 
spread observed in the old stellar population of 
dSphs is not a definite signature of star formation 
extended over 2-4 Gyrs (Ikuta & Arimoto 2002; 
Grebel & Gallagher 2004). Neither the argument 
based on the long timescale {e.g., 1-4 Gyr) needed 
to have a substantial iron production by Type la 
SN is a strong argument in favor of prolonged star 
formation. The typical time scale for iron enrich- 
ment by Type la SN is uncertain and very likely 
is dependent on the mode of star formation. Mat- 
teucci & Recchi (2001) estimated that this time 
scale varies from 40-50 Myr for an instantaneous 
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Fig. 10. — Metallicity spread of the stellar popu- 
lations <J[Fe/H] 1 8.S a function of mean metallicity 
[Fe/H] for the simulated (gray squares) and ob- 
served dSphs (symbol markings are the same as 
in Fig. 6). The spread (J[pe/H] is defined as the 
variance of the number weighted metallicity dis- 
tribution of the "stellar particles" in each dwarf. 
A rough estimate of the uncertainty on the metal- 
licity spread is shown by the circled errorbar on 
the top left corner. The error is estimated as the 
mean difference between (J[Fe/H] using the mass 
weighed and the number weighed metallicity dis- 
tributions. 

starburst to 0.3 Gyr for a typical elliptical galaxy 
to 4-5 Gyr for a disk of a spiral galaxy like the 
Milky Way. The often quoted 1 Gyr time scale is 
based on observations of [a/Fe] versus [Fe/H] in 
the solar neighborhood stars and is not an univer- 
sal value. A more detailed study would be needed 
to derive accurate metal abundance ratios based 
on the star formation and merger histories found 
in the simulations. This is beyond the aim of the 
present paper but it is worth pursuing in a sepa- 
rate study. 

We can also study the stellar metallicity distri- 
bution in each one of our simulated dwarfs. Fig- 
ure 11 show such an example. The solid and 
dashed lines show the number and mass weighted 
distributions respectively. The metallicity distri- 




[Fe/H] 

Fig. 11. — The stellar metallicity distribution of 
one simulated galaxy. The solid and dashed lines 
show the number and mass weighted distributions 
respectively. 

bution is skewed toward low metallicities but only 
a small fraction of the stars, having metallicities 
< 10~^ Zq, could be defined as Population III 
stars. 

4. Discussion 

In this paper we have proposed that dwarf 

galaxies of the Local Group (and, by implication, 
all other dwarf galaxies in the universe) formed in 
three different evolutionary paths: 

"true fossils" formed most of their stars in the 
prc-rcionization era, and have had little (say, 
less than 30%) star formation since then; 

"polluted fossils" started as true fossils, but 
had substantial episode(s) of subsequent star 
formation as they continued accreting mass 
and were tidally shocked during the forma- 
tion of the parent galaxy halo; finally 

"survivors" started forming stars mostly after 
reionization. 

This conclusion is motivated by considering star 
formation histories of Local Group dwarfs and by 
comparing their properties with properties of sim- 
ulated galaxies that formed all their stars before 
cosmological reionization. The fact that our simu- 
lated galaxies are remarkably similar to the subset 
of the Local Groups dwarfs that we identified as 
"true fossils" in Figure 5 renders support to our 
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conclusions. This subset includes about all known 
dSphs. 

To avoid confusion it is worth noting that the 
structural properties of the galaxies that we call 
"true fossils" are the result of feedback processes 
in action prior to reionization. The effect of reion- 
ization is simply to preserve those properties by 
suppressing the star formation rate in halos with 
masses that remain smaller than the filtering mass 
in the IGM. In this sense our model differs from 
the standard "reionization" scenario for the for- 
mation of dwarf galaxies. 

How would one test this hypothesis? For exam- 
ple, the main argument of Kravtsov ct al. (2004) 
for the tidal origin of the majority of Local Group 
dwarfs is their preferential location close to the 
parent galaxy. However, true fossils, living in the 
oldest dark matter halos of their mass, are highly 
biased, and are also preferentially located closer to 
the center of the parent halo. Moreover there are 
at least two examples of dSph galaxies that are 
quite far from massive galaxies: Cetus and Tu- 
cana, and perhaps the newly discovered Applesl 
(Pasquah et al. 2005). 

A more powerful test is based on our conclusion 
that "survivors" form later than the "true fossils" . 
While measuring absolute ages of dwarf galaxies to 
1 Gyr precision is not possible at the moment, de- 
termining the relative difference of about 1 Gyr in 
age between two galaxies observed by a uniform 
technique may, in fact, be possible for a specific 
range of mctallicitics and ages using the so-called 
"horizontal branch index" (G. Gilmore, private 
communication; also see Zinn 1993; Harbeck et al. 
2001; Makcy & Gilmore 2004): younger galaxies 
have redder horizontal branch for the same metal- 
licity, although, how precise the horizontal branch 
index method is and the role of the so-called "sec- 
ond parameter effect" is still a matter of debate 
(the anonymous referee, private communication). 

Luckily, some of the Local Group dwarfs fall 
into that range of parameters. For example, two 
galaxies of our "true fossil" group - And I and And 
H - have [Fe/H] « —1.5, similar to the iron abun- 
dance in two of "survivors" - WLM and NGC3109 
(having pairs of galaxies with the same metallic- 
ity is absolutely crucial as even small differences in 
metallicities can obscure large differences in age). 
Thus, we can make a prediction that most of the 
stars that formed in the first burst of star for- 



mation in WLM and NGC3109 should be about 
1.5 Gyr younger than most of stars in And I and 
And II. The color-magnitude diagrams of And I, 
And II, and WLM are measured down toV> 26, 
well below the horizontal branch (Da Costa, Ar- 
mandroff, & Caldwell 2002; Rejkuba et al. 2000), 
although, to the best of our knowledge, no pho- 
tometry down to that levels exists for NGC3109. 
Visual comparison of the horizontal branch mor- 
phologies for the first three galaxies indicate that 
And I and And II are indeed substantially older 
than WLM, also a more rigorous, quantitative test 
is required to confirm the visual impression. Such 
a comparison would be a crucial test of our hy- 
pothesis - if it is found that And I and And II are 
not substantially older than WLM and NGC3109, 
that would imply that they are not "true fossils" , 
and by inference, few other dwarf spheroidals are. 
The main problem that affects this method when 
comparing the old populations of dSphs to dirrs is 
that normally dIrrs have mixed populations and 
the overlap between the blue horizontal branch 
and the young main sequence make it difficult to 
derive reliable horizontal branch morphology in- 
dices in dIrrs. Perhaps a more reliable and ac- 
curate relative age indicator is the main sequence 
turn-off. The internal accuracy can be better than 
1 Gyr in few nearby dwarfs with photometric data 
that extends below the main-sequence turnoff. 

Using available data, Grebel & Gallagher 
(2004) address the impact of reionization on the 
star formation histories of dwarf galaxies. They 
conclude that most or all dwarf galaxies present an 
old population that is coeval with the oldest Milky 
Way populations. They also find that all of the 
Local Group dwarf galaxies exhibit considerable 
abundance spreads even within their old popula- 
tions indicating that star formation did not occur 
in a single burst, but in episodes that must have 
extended over several Gyr, well beyond reioniza- 
tion. In this paper we have argued that the ob- 
served metallicity spreads are not necessarily an 
indication of star formation extended over several 
Gyr. In our simulations all the dwarfs that form 
before reionization have metallicity spreads consis- 
tent with observations. The large iron abundances 
in dSphs have also been used as an argument in 
favor of an extended star formation. This is be- 
cause the time-scale for iron enrichment is often 
assumed to be longer than 1 Gyr. But we point out 
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that this assumption may be true for a continuous 
mode of star formation, but is not necessarily uni- 
versally true. The time scale depends on the star 
formation history and could be as short as 50 Myr 
for a bursting mode of star formation (Matteucci 
& Recchi 2001). 

In principle, [a/Fe] ratio could be used to dis- 
criminate between star formation histories consist- 
ing of a single burst or a more prolonged episode 
of star formation. At the moment it is not clear 
how much this test could distinguish "true fossils" 
from "survivors" , since both are expected to have 
the main burst of star formation for a relatively 
short period of time. The star formation rate in 
the "true fossils" should decrease substantially af- 
ter reionization. In addition our simulations show 
that they have a very peculiar mode of star forma- 
tion characterized by multiple short bursts (about 
10 — 100 Myr long). It is unclear to us how this 
mode of star formation may affect their [a/Fe] ra- 
tio and be used as a discriminant. Using exist- 
ing data on [a/Fc] ratios in dSphs (Tolstoy et al. 
2003, Venn et al. 2004) it might be possible to 
distinguish between different star formation and 
enrichment histories but this would require a more 
detailed modelling that is beyond the aim of the 
present paper. 

In addition, we find that the most massive end 
of the "true fossils" population in our simulation 
(with mass M > 1 — 2 x 10* Mq) retain a fraction 
of their gas even after reionization (see Fig. 1). 
These few "transition" objects may continue form- 
ing stars for a while after reionization and have a 
more extended period of star formation with re- 
spect to the rest of the "fossil" population. 

To conclude, we would like to emphasize an- 
other exciting possibility. From the results of our 
simulations, we feel quite confident in predicting 
the existence of fainter and lower surface bright- 
ness galaxies than the one currently discovered 
(see the gray squares at the extreme left corners 
in Figs. 3, 6 and 7). Two promising candidates for 
such a class of dwarf galaxies were recently discov- 
ered in M31 and near the Milky Way (Zucker et 
al. 2004, Willman et al. 2005). 

As we mentioned above, it is mostly likely that 
the majority of reionization fossils will merge into 
larger objects by z = 0, and only a small fraction 
of them will survive. Thus, we can introduce a 
"survival probability" S (which can be a function 
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Fig. 12.— Prediction of the number of faint dSphs 
per logarithmic luminosity bin in the Local Group 
(histograms). The points with error bars (Poisson 
errors) show the number of observed dwarf galax- 
ies in the Local Group. The solid line refers to the 
case in which the survival probability of reioniza- 
tion fossils is constant as a function of luminosity, 
while the dashed line shows an example of a case in 
which the survival probability is mildly luminosity 
dependent (see text for details). 

of other parameters, in particular the luminosity) 
of a reionization fossil to remain intact by the 
present time. This probability can be estimated by 
comparing the number of observed dwarfs in the 
Local Group (the total mass is about 4 x 10^^ M©) 
with the number of simulated galaxies in our vol- 
ume (with the total mass of 1.2 x 10^^ M©). If 
we assume that both the observations and our 
simulation are complete for galaxies with Ly be- 
tween 10^ Lq and 3 x 10^ L©, we find a value of 
S = 3% for the survival probability in this lumi- 
nosity range. 

Figure 12 shows our predictions for the abun- 
dance of these "ultra-dwarfs" in the Local Group. 
The solid line in Fig. 12 refers to the case in 
which the survival probability of reionization 
fossils is kept constant at 3% - in this case; 
the observations may be missing about half 
of the satellites with luminosities in the range 
2 X lO-"^ Lq < iv- ^ 6 X lO-"^ L© (due to a small 
size of the computational box, out simulation 
is incomplete for Ly ~ 10^ L© , as we discussed 
above). Alternatively, we can assume that the ob- 
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servations are complete in this luminosity range, 
in which case the survival probability must be a 
weak function of luminosity. Fitting the power- 
law dependence (as an example) to the abundance 
of observed dwarfs in the three lowest luminos- 
ity bins (2 X 10^ Lq < < 6 X 10^ L©), we find 
S{Lv) = 3% X {Lv/2 X 10^ Lq)^/^ (this fit reaches 
S = 100% for the Milky Way luminosity, and so, 
at least, makes sense). This case is shown with the 
dashed line in Fig. 12. As one can see, a modest 
change in the survival probability makes a large 
change in the predicted numbers of dwarf galaxies 
in the Ly ^ IQ^ luminosity range, and, thus, 
predicting the precise abundance of "ultra-dwarfs" 
will be difficult. On the other hand, measuring 
the abundance of dwarf galaxies in this luminos- 
ity range will put non-trivial constraints on the 
luminosity dependence of the survival probabil- 
ity, and, therefore, on the hierarchical merging 
process and galaxy clustering on extremely small 
mass scales. 

Although the detection of such objects in the 
Milky Way is very difficult if not impossible us- 
ing photometric techniques (due to contamination 
from foreground stars), in the future it might be- 
come possible to search for these objects using 
spectroscopic techniques by measuring the radial 
velocities of many stars. Any observation able to 
put constraints on properties of the faintest bound 
object in the Local Group would be of extreme in- 
terests not only for understanding galaxy forma- 
tion in the early universe, but, more importantly, 
for uncovering the properties of the elusive dark 
matter particles. 
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